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Abstract 

When it is oxidized, graphite can be easily exfohated forming graphene oxide (GO). GO is a 
critical intermediate for massive production of graphene, and it is also an important material with 
various application potentials. With many different oxidation species randomly distributed on 
the basal plane, GO has a complicated nonstoichiometric atomic structure that is still not well 
understood in spite of of intensive studies involving many experimental techniques. Controversies 
often exist in experimental data interpretation. We report here a first principles study on binding 
energy of carbon Is orbital in GO. The calculated results can be well used to interpret experimental 
X-ray photoelectron spectroscopy (XPS) data and provide a unified spectral assignment. Based on 
the first principles understanding of XPS, a GO structure model containing new oxidation species 
epoxy pair and epoxy-hydroxy pair is proposed. Our results demonstrate that first principles 
computational spectroscopy provides a powerful means to investigate GO structure. 
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I. INTRODUCTION 



Graphene is one of the most studied materials in recent years owing to its peculiar physical 
properties and great potential for various applications. A big technological challenge 

in graphene research is to massively produce high quality samples. Epitaxial growth rep- 



resents a suitable way for electronics application. [4, [s], Q] However, for general purpose 
applications, solution based strategy is more attractive. One promising graphene synthesis 
route is oxidizing and exfoliating graphite, then reducing the obtained graphene oxide (GO) 



sheets. 



fly, 



Id, 
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12| By this way, the pristine graphene properties can be largely restored 



during the reduction. At the same time, GO itself is an important material with various 



applications. It is a good candidate for 
shows excellent mechanic properties. 



composite materials. 



13| and its hydrate, GO paper. 



14| It is possible to integrate GO into all-graphene 



electronic devices [15|] and to obtain near-UV to blue photoluminescence from GO. ^] GO 



based hydrogen storage materials are suggested to be very promising. [17|, [181, [lOj By us- 
ing GO support, Scheuermann et al. [20] also achieved much higher Pd catalysis activity 
compared to conventional Pd/C system. 

Because of its importance, it is highly desirable to understand the atomic structure of 
GO and the oxidation processes to prepare it. There is a long history of GO structure 



studies. 



21 



22 
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26| In the most popular structure model, GO is constituted 



of epoxy and hydroxy groups attached to the carbon layer, and terminated with hydroxy 

and carboxyl groups. Many experimental techniques, such as nuclear ma gnet ic resonance 

(NMR), [25, 



26 



27 



28| x-ray photoemission spectroscopy (XPS), 



26 



29 



vibrational spectroscopy. 



33, 



31 



3^ and 



3J] have been employed to explore the structure information. 



However, the interpretation of experimental data turns out to be difficult, since GO is a 
nonstoichiometric compound with randomly distributed oxidation groups. 

XPS is one of the best local structure probes to identify oxidation species on GO. Ex- 
perimental XPS for heavily oxidized GO have mainly three peaks, namely Pi, P2, and P3 
from lower to higher binding energies. It is commonly accepted that Pi comes from sp'^ car- 
bon, and P2 is contributed by carbon atoms connected to epoxy (C-E) and hydroxy (C-OH) 
groups. However, definitive assignment for P3 is currentl^Mtiot available. Szabo et al. 26] 
assigned this peak to ketones (C=0), while Jeong et al. [29] attributed it to edge carboxyl 
(COOH) groups. Considering its high intensity, peak P3 can not be solely contributed by 
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edge groups. On the other hand, the existence of interior ketone species has been excluded 
by other experiments. For example, a recent two dimensional NMR experiment suggested 
;hat carbonyl groups are spatially separated from the majority sp^, C-OH, and C-E carbons. 



27| This result implies that C=0 groups mainly distribute at GO edges, as also suggested 



by a recent O Is X-ray absorption spectrum measurement. 35| Therefore, both experimen- 
tal XPS assignments are questionable, and new species may exist in GO [S] that has a 
contribution to P3. 

Recently, new functional groups such as epoxy pair (EP) and manganate ester (C-MnO^ 
have been proposed in studies on oxidation induced cut of graphitic materials. 36|, l37 | 



When epoxy chain breaks underlying C-C bonds, [38| additional epoxy groups can attach 
to it forming EP groups. EP can then transit to carbonyl pair (CP) and finally break the 
carbon network. However, there is a significant dissociation barrier for the first EP to CP 
transition, 36|] which makes EP can in principle be widely existed in GO. Besides these two 



groups, to reach a tentative assignment of NMR spectrum, a new structure motif of lactol 
rings is suggested to be existed at GO edge sites. (28| However, XPS signatures of these 
new spices are still unclear. 

Theoretical simulation can provide useful information to obtain correct spectroscopic as- 
signments and thus a more reliable GO structure model. Unfortunately, previous theoretical 



studies on GO mainly focus on energetics 



36 



39, 



40 



41 



42l | instead of properties that can 



be directly compared with experimental data. The complexity of the GO potential energy 
surface limits the power of energetics study to some extent, especially when artificial pe- 
riodic boundary condition must be adopted. Kudin et al. 33|] have simulated the Raman 
spectra of GO, and proposed an alternating single-double carbon bond model. Compared 
to Raman spectroscopy, XPS is a more localized probe and thus a better choice to study 
local structures. 

In this study, we calculate core chemical shift of carbon Is orbital for possible oxidation 
species in GO using density functional theory, which can be compared to experimental 
XPS spectra directly and leads to a unified reassignment. Based on such a computational 
spectroscopy strategy, we propose a GO structure model with new species epoxy pair and 
hydroxy-epoxy pair. The remainder of this article is organized as follows. In section 2, 
we introduce the computational methods. XPS simulation results and implications to GO 
structure are present in section 3. Finally, we conclude in section 4. 
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II. COMPUTATIONAL DETAILS 



XPS simulation with plane wave basis set. Spin polarized electronic structure 
calculations are performed with the Vienna ab-initio simulation package (VASP) 43|] using 
Perdew, Burke and Ernzerhof (PBE) 4J] exchange-correlation functional. The electron-ion 
interaction is described with projector augmented wave (PAW) method [45], which permits 
to set the occupation of core level as any number between and 1. Different two dimensional 
sheet and one dimensional nanoribbon models are adopted to simulate GO. The binding 
energy of C Is orbital is calculated as the energy difference between the ground state and 
the core-excited state with one core electron removed. 46| We have tested this computational 
approach by calculating the Is core orbital binding energies of C, N, and O atoms in alanine 
and tryptophan, and the obtained core chemical shifts in different chemical environments 
agree well with the corresponding experimental values. 47] 

XPS simulation with atomic basis set. Cluster models are also used to calculated 
he binding energy of C Is orbital. The geometries of all clusters are optimized with DMoP 



49| and the binding energy calculations are performed with StoBe-deMon. [50|] In 



binding energy calculations, gradient corrected Becke (B88)[51| exchange functional and 
Perdew (P86)|52| correlation functional are adopted. Effective core potential is used to 
describe the non excited atoms and for ionized atoms, lGLO-111 basis set [53] is used instead. 
To get the binding energy of core electron, one C Is electron is removed from the system. 
The core chemical shift is defined as the energy difference between core electron removed 
system and the ground state. Such a XPS simulation protocol has been widely used in 



54 



55 



56| For 



many studies, and it generally leads to a good agreement with experiment, 
different cluster models, the core chemical shift is calibrated according to the binding energy 
of C-H far way from the oxidation groups. 



III. RESULTS AND DISCUSSION 



The main task of this study is to calculate the core chemical shift of C Is orbital for 
different oxidation species. For this purpose, two first principles simulation protocols and 
many different structure models are used. The absolute value of core chemical shift can be 
affected by many factors, for examples, by charging effect in experiment and by work function 
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of the theoretical GO model. We find that for a perfect graphene sheet, the work function 
has a value of 4.7 eV, while for a hydrogen saturated ID graphene nanoribbon, it becomes 
3.7 eV. It indicates that the structure model used for GO can have a significant effect on 
the work function and thus on the absolute value of core chemical shift. Therefore, we only 
focus on relative core chemical shift (R-CCS). The average binding energy of epoxide group 
in each structure model is set as a common reference to compare results from calculations 
and from experiments. 




FIG. 1: Cluster models used to calculate C Is binding energies, (a) epoxy, (b) two OH groups, 
(c) two edge C=0 groups, (d) carboxyl group, (e) line defect epoxy group (C-E') and epoxy pair 
(C-EP) group, (f) C-E' and carbonyl pair (CP), (g) five-ring lactol group, (h) six-ring lactol group, 
(i) hydroxy-epoxy pair, and (j) Mn04- attached structure. 
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We first consider the well studied epoxy and hydroxyl groups, which contribute to XPS 
peak P2. In some experiments, P2 is further deconvoluted into two sub-peaks, namely 



and P|. 



29 



3l| The former is about 0.8 eV lower in binding energy than the latter. These 
two peaks have been tentatively assigned to C-OH and C-E, respectively. However, such 
an assumption needs more stringent test, for example, by first-principles calculation on the 
binding energy of carbon Is orbital. 
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FIG. 2: (a) Graphene sheet with isolated epoxy and hydroxy groups. Dashed lines mark a (6 x 6) 
unit cell, (b) A GO model with epoxy and hydroxy groups. Dashed lines indicate the (4 x 4) cell 
used in our calculations. Red and green spheres represent oxygen and hydrogen atoms, respectively. 



As shown in Figure [T^ and[T]D, we construct two clusters, one with C-E and the other with 
C-OH. Both clusters have an exclusive Clar's structure, where all vr-electrons form separated 
aromatic sextets. 57| Dangling bonds are saturated with hydrogen, and the oxidation groups 
are put in the center of the cluster. The calculated R-CCS of C-OH with respect to C-E 
is only 0.1 eV, indicating that C-E and C-OH are difficult to be resolved in XPS. It is also 
interesting to note that the core binding energy of C-OH is larger than that of C-E, contrary 
to the tentative assignment for P2 and P2. 

For periodic systems, we put two isolated hydroxy and epoxy groups in a big (6 x 6) 
supercell, as shown in Figure [2^. R-CCS of sp^ C with respect to C-E is -1.7 eV, while the 
difference between the binding energies of C-E and C-OH is still within 0.1 eV, consistent 
with the atomic basis set calculation. To take the coupling between epoxy and hydroxy 



groups into account, we also consider a more realistic GO structure model 



40j (Figure [2)d). 



In this case, the binding energy of C-OH is larger than that of C-E by 0.3 eV, showing 
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the important effect of chemical environment. However, such a small difference is still not 



distinguishable in XPS experiments, 
required. 



SOj Therefore, a new assignment for P2 and P| is 



' 4 \ 

(a) (b) 



FIG. 3: New oxidation species, (a) epxoy pair, (b) hydroxy-epxoy pair. The corresponding car- 
bon atoms (red) are named as C-EP and E-C-OH, respectively. Both C-EP and E-C-OH are 
monovalently connected to two oxygen atoms. 

To better understand the XPS peaks, more oxidation species beyond C-OH and C-E 
should be considered. First, we take edge groups. It is known that there are C=0, COOH, 
and hydroxy (C-OH-edge) groups existing at GO edge. In this study, we also consider other 
edge groups including C-H, CH3, CO3, and lactol (Figure [1^ and [T)i). For new interior 



groups, some clues can be found from the oxidation induced cutting of graphene. [36| In 



highly oxidized GO, some epoxy groups may align in a line, which leads to the break of the 



underlying C-C bonds. 38|] Then, a new epoxy group can be attached to form a EP group 
(C-EP, Figure ISi), 3^ or a hydroxy can be added to from a epoxy-hydroxy pair (E-C-OH, 



Figure Hb). E-C-OH is actually also existed in the recently proposed edge lactol groups, [28 1 
and the latter can be considered as a special case of the former. 

Once the wide existence of aligned epoxy groups is recognized, the formation of C-EP 
and E-C-OH becomes very natural. When aligned epoxy groups break the underlying C-C 
bonds, the corresponding carbon atoms {a carbon) are connected to two other carbon atoms 
{/3 carbon) and the epoxy group, to which the a carbon is bound with one valence electron. 
Therefore, there is a tt electron left on the a carbon. If both the two (3 carbon atoms are 
oxidized to be sp^ hybridized, then the a carbon becomes a radical. In highly oxidized 
samples, the possibility that (3 carbon atoms are saturated is quite high. Therefore, the a 
carbon is very possible to be attacked to form epoxy pair or epxoy-hydroxy pair groups in 
the oxidizing environment. 

Together with C-E and C-OH, these new oxidation species have been put in various of 
structure models. Individual groups have been put on clusters (Figured]). We also construct 



FIG. 4: (a) A GO nanoribbon model with zigzag edges. The part between two dashed Hnes is a 
unit cell along the ribbon direction, (b) A GO structure model based on epoxy line defect. 



several GO nanoribbon models with many oxidation species randomly distributed both on 
the plane and at the edges. One of these ribbon models with zigzag edges is shown in 
Figure Hti. It turns out that the calculated core chemical shifts are almost not dependent 
on whether zigzag or armchair edges are adopted in the GO model. Another kind of GO 
models we considered is epoxy line defect based, for which different functional groups are 
put on the epoxy chain, as shown in Figure |Dd. 

With these GO models, we have calculated the core chemical shift for oxidation species C- 
OH, C-E, C-EP, and E-C-OH on planes, and C-OH-edge, C=0, C-H, CHg, COOH, COOO, 
and lactol at edges. The calculated R-CCS values are collected in Table [U together with 



experimental values reported in the literature. 



26 



291] Chemical shifts for all calculated 
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oxidation species are well within the energy window of the experimental XPS data. Again, 
we find that epoxy and hydroxy groups in the interior of ribbons are not distinguishable. 
The peak P2 originally assigned to C-OH by Jeong et al. 29| has contributions from edge 
attached C-OH-edge and C=0. 

We notice that, in experimental spectra, have comparable intensity to P^- Therefore, 
there should be other contributions to P2 besides the edge groups. For epoxy chains with 
the underlying C-C bonds broken, the corresponding carbon atoms (C-E') almost keep sp^ 
hybridization with a C-0 bond length about 1.37 A, very close to that in C-OH-edge (1.36 
A). R-CCS of C-E' is also very similar to that of C-OH-edge. Since it is distributed in the 
interior part of GO, C-E' can make a strong peak in XPS. 

Our calculations show both C-EP and E-C-OH contribute to the XPS peak P3. R-CCS 
of C-EP in an epoxy line defect (1.5 eV) is slightly higher than its isolated counterparts (0.9 
to 1.0 eV). E-C-OH has a larger core level binding energy compared to C-EP. Edge group 
CO OH has a similar R-CCS with C-EP. Since the newly proposed five- or six-membered 
lactol ring contain E-C-OH, its R-CCS also locates in the region of P3, only slightly lower 
than that of interior hydroxy-epoxy pair. However, as an edge groups like COOH, it only has 
limited contribution to XPS. Main contribution to P3 thus comes from C-EP and E-C-OH. 

EP can dissociate to form a carbonyl pair (CP, Figure [T]F andHb). 36(] R-CCS of an 
individual CP is close to zero, significantly higher than edge C=0 group, indicating its 
chemical environment is strongly affected by neighboring C-E'. If the whole line defect is 
broken, CP becomes C=0 edge groups with a negative R-CCS. Since an existing CP strongly 
lowers the dissociation barrier of neighboring EPs, isolated CPs are not expected to be widely 



existed. 



361 ] However, an important fact disclosed here it that, even if C=0 groups exist in 



;he interior of GO, its contribution to XPS goes to P2 instead of the originally supposed P3. 



26| 



Another spices mentioned in the literature is C-MnOj. In our VASP calculations with 
plane wave basis set, for both a (6x6) unitcell and a cluster model, we fail to get a local 
minimum structure with Mn04 adsorbed on graphene, MnOj always moves away during 
the optimization. This result suggests that C-Mn04 may be not stable or only stable in 
solution, and thus does not exist during XPS measurements. An adsorbed structure can 
be obtained using the atomic basis set with the DMoF package. With this geometry, the 
calculated R-CCS is about -0.7 eV, far away from peak P3. 
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TABLE I: R-CCS of carbon Is orbital for various oxidation species reported in experiment and 
calculated with atomic and plane wave basis sets in this work. The binding energy of C-E (287.2 eV) 
is set as a common reference point. The numbers in parentheses are the linewidths of experimental 
peaks. 





Atomic 


Plane Wave 


EXP[26J 


EXP [29] 


C-H 


-1.8 


-2.5 ~ -1.9 






C sp'^ 


-1.7 


-2.3 ^ -1.6 


-3.1 


-2.2 (±1.2) 


CH3 




-1.8 ~ -1.6 






C-OH-edge 




-0.9 ~ -0.6 




-0.8 (±1.0)'= 


C=0 


-0.7 


-0.9 ~ -0.1 




-0.8 (±1.0)^ 


C-Mn04 


-0.7 








C-E' 


-0.7~-0.5 


-0.8 ^ -0.3 






C-OH 


+0.1 


-0.1 ^ +0.3 


0.0 (±1.2) 




C-E 


0.0 


-0.3 ^ +0.3 


0.0 (±1.2) 


0.0 (±1.0) 


CP 


0.0~+0.1 


-0.3 ~ -0.2 






C-EP 


+1.1 


+0.9 ~ +1.5 


+2.1 (±1.2)" 


+1.8 (±1.2)^^ 


COO(H) 


+1.3~+1.9 


+0.9 ~ +1.4 


+2.1 (±1.2)" 


+1.8 (±1.2)'^ 


E-C-OH 


+1.7~+1.9 


+1.5 ~ +1.9 


+2.1 (±1.2)" 


+1.8 (±1.2)'^ 


COOO 




+3.0 


+3.6 (±1.0)^ 



" assigned to carbonyl groups in Ref. |26]. 
^ assigned to carboxyl groups in Ref. 
assigned to hydroxy groups in Ref. 
assigned to carboxyl groups in Ref. [29] 



In some experimental XPS data, there is a small peak with a R-CCS more than 3 eV. 26| 
Our calculations suggest that it can be related to carbon atom connected to three oxygen, 
possibly at GO edges. Until now, as shown in Figure [5l we reach a unified assignment 
for experimental XPS. The lowest peak Pi comes from sp^ carbon. Peak P2 has two con- 
tributions: the lower one comes from edge groups and edge-like group C-E', and the higher 
one is contributed by epoxy and hydroxy groups. Peak P3 mainly comes from C-EP and 
E-C-OH, with also small contributions from edge COOH groups. 
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FIG. 5: Schematic experimental XPS curve and our first principles assignment. 
With the proper first-principles assignments, the XPS intensity evolution observed in 



experiments 



26 



29| can also be well understood. First of all, when the degree of oxidation 



increases, the peak of sp^ carbon (Pi) should decrease or even disappear for completely 



oxidized samp 
per iment ally. 



es. At the same time, an increase of P3 compared to P2 was observed ex- 



26 



291] which is because more EP and E-C-OH groups appear during the 
I ! 

oxidation. Jeong et al. ^] observed an increase of P2 and a decrease of P2 during oxida- 
tion, which can also be understood with our assignments. In fact, during the increase of 
the degree of oxidation, the possibility for epoxy groups to align in a line increases, which 
leads to more C-E' groups (^2")- When GO is heated, a significant decrease of P3 has been 
observed. 29|] The dissociation of EP to CP should have made an important contribution 
to it. 



IV. CONCLUSIONS 

Our DFT calculations provide relative C Is core chemical shifts for possible oxidation 
species on GO, whose binding energy follows the order: C sp^ < edge groups (such as C=0, 
C-OH-edge, and C-E') < C-OH and C-E < C-EP, COOH, and E-C-OH < COOO. These 
results lead to a reassignment of experimental XPS spectra: (1) P2 is mainly contributed 
by C-E and C-OH, and they are not distinguishable; (2) edge groups and edge-like C-E' 
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groups has lower binding energies compared to other P2 groups; (3) although it has been 
previously used to explain P3 in XPS experiments, C=0 actually has a low binding energy 
corresponding to P2; (4) -P3 mainly comes from new species C-EP and E-C-OH; (5) even 
higher binding energy may obtained from COOO. 

Based on our first principles spectroscopic understanding, we reach the following two 
important conclusions about the atomic structure of GO: (1) epoxy and hydroxy groups can 
be closely packed together forming epoxy pairs and epoxy-hydroxy pairs in highly oxidized 
samples; (2) as the main species contributed to the lower energy part of P2, aligned epoxy 
group (C-E') represents an important GO structure motif, which should be considered in 
all GO property studies and not only in oxidation induced cutting. Theoretical simulations 
also lead to an understanding of the experimental evolution of XPS for different samples, 
which provides us some insights to the oxidation processes of GO. Our results demonstrate 
that first principles computational spectroscopy is a very promising tool in nano structure 
research. 
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